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Abstract. - We perform comprehensive theoretical analysis of high magnetic field behavior of the 
heavy-fermion (HF) compound YbRli2Si2. At low magnetic fields B, YbRli2Si2 has a quantum 
critical point related to the suppression of antiferromagnetic ordering at a critical magnetic field 
B _L c of B = B c o — 0.06 T. Our calculations of the thermodynamic properties of YbRli2Si2 
in wide magnetic field range from B c o — 0.06 T to B ~ 18 T allow us to straddle a possible 
metamagnetic transition region and probe the properties of both low-field HF liquid and high-field 
fully polarized one. Namely, high magnetic fields B ~ B* ~ 10 T fully polarize corresponding 
quasiparticle band generating Landau Fermi liquid (LFL) state and suppressing HF (actually 
NFL) one, while at elevating temperatures both HF state and corresponding NFL properties 
are restored. Our calculations are in good agreement with experimental facts and show that 
the fermion condensation quantum phase transition is indeed responsible for the observed NFL 
behavior and quasiparticles survive both high temperatures and high magnetic fields. 
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An explanation of the rich and striking behavior of 
heavy fermion (HF) metals is, as years before, among the 
main problems of modern condensed matter physics. One 
of the most interesting and puzzling issues in the research 
of HF compounds is their non- Fermi liquid (NFL) behav- 
ior in a wide range of temperatures T and magnetic fields 
B. For example, recent measurements of the specific heat 
C of YbRli2Si2 under the application of magnetic field B 
show that the above temperature range extends at least 
up to twenty Kelvins as reported in the inset to fig. [U 
As it is well-known from Landau Fermi liquid (LFL) the- 
ory, the ratio C/T is proportional to quasiparticle effective 
mass M*. The inset to fig. Q] reports the dependence of 
C(T)/T, which has a maximum M^ 13X (B) at some temper- 
ature T max (B). It is seen from the inset, that M* lax (B) 
decreases as magnetic field B grows, while T max (T?) shifts 
to higher T reaching 15KatB = 18T[T]. 
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A deeper insight into the behavior of C/T in the inset 
to fig. Q] can be achieved using some "internal" scales. 
Namely, near QCP it is convenient to divide the effec- 
tive mass M* and temperature T by their maximal val- 
ues, M,* ax and T max respectively. This generates the nor- 
malized effective mass = M*/M^ ax and temperature 
Tjy = T/T max [2]. In the main panel of fig. Q]the obtained 
dependence Mjy(Tjv) is shown by symbols, corresponding 
to different magnetic fields. This immediately reveals the 
scaling in the normalized experimental curves - the curves 
at different magnetic fields B merge into a single one in 
terms of the normalized variable TV = T '/T max . It is seen 
from fig. \T\ that the normalized effective mass M^(TV) 
is not a constant as it would be for LFL case. Rather, it 
shows the scaling behavior in normalized temperature T/y. 
It is also seen from fig. [1] (both the main panel and inset) 
that the NFL behavior and the associated scaling extend 
at least to temperatures up to twenty Kelvins. 

Thus, we conclude that a challenging problem for the- 
ories considering the high magnetic field ( B ~ B*) NFL 
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behavior of the HF metals is to explain both the scaling 
and the shape of M^(T/v). Another part of the problem 
is the remarkably large temperature and magnetic field 
ranges where the NFL behavior and scaling are observed. 
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Fig. 1: The normalized effective mass = M* /M^ax versus 
normalized temperature Tm — T/T max . is extracted from 
the measurements (shown in the inset) of the specific heat C/T 
on YbRh2Si2 in magnetic fields B Q] listed in the legend. Our 
calculations (made at B ~ B* when the quasiparticle band 
is fully polarized) are depicted by the solid curve tracing the 
scaling behavior of Mjy. The inset reports the temperature 
dependence of the electronic specific heat C/T of YbRli2Si2 at 
different magnetic fields [T] shown in the main panel legend. 
The illustrative values of Af,* ax and T max at B — 8 T are also 
shown. 



In this letter, based on the theory of fermion conden- 
sation quantum phase transition (FCQPT) [5] we analyze 
the thermodynamic properties of YbRh.2Si2 at both low 
and high magnetic fields. Our calculations of the specific 
heat and magnetization allow us to conclude that under 
the application of magnetic field the heavy-electron system 
of YbRIi2Si2 evolves continuously without a metamagnetic 
transition. At low temperatures and high magnetic fields 
B ~ B* the system is completely polarized and demon- 
strates the LFL behavior, while at elevated temperatures 
the HF behavior and related NFL one are restored. The 
obtained results are in good agreement with experimen- 
tal facts in the entire magnetic field (0.1 T - 18 T) and 
temperature (40 mK - 20 K) domains. 

In our FCQPT approach [2J, to study the (generally 
speaking NFL) behavior of the effective mass M*(T,B), 
we simply use Landau equation for the quasiparticle effec- 
tive mass in a Fermi liquid. The only modification is that 
in our formalism the effective mass is no more constant but 
depends on temperature, magnetic field and other exter- 
nal parameters. For the model of homogeneous HF liquid 
at finite temperatures and magnetic fields, this equation 



acquires the form J3J[3] 

1 = 1 + E 

rn ^ — * 



M*{T,B) 



P£R F 

P F 



(Pf,p) 



dn ai {p,T,B) dp 



dp 



(2tt)= 



(1) 



where m is a bare electron mass, F - jCri (pp,p) is the Lan- 
dau amplitude, which depends on Fermi momentum pp, 
momentum p and spin a. Here we use the units where 
H = ks = 1- For definiteness, we assume that the HF 
liquid is 3D liquid. The Landau amplitude has the form 
® 
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where E[n] is the system energy, which is a functional of 
the quasiparticle distribution function n [2J[3]. It can be 

expressed as 



n a (p,T) = < 1 + cxp 
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where e(p,T) is the single-particle spectrum. In our case, 
the chemical potential p, depends on the spin due to Zee- 
man splitting p, a = p ± p-bB, p,B is Bohr magneton. 

In LFL theory, the single-particle spectrum is a vari- 
ational derivative of the system energy E[n a {p 1 T)\ with 
respect to occupation number n, e(p,T) = SE[n(jp)]/Sn. 
Choice of the amplitude is dictated by the fact that the 
system has to be at the quantum critical point (QCP) of 
FCQPT. Namely, in this region the momentum-dependent 
part of Landau amplitude can be taken in the form of 
truncated power series F — a(p — p') 2 + b(p — p') 3 + 
c(p — p') 4 + where a, b and c are fitting parameters. 
We note that this interaction, being an analytical func- 
tion of (p — p') 2 , can generate topological phase transi- 
tions interfering in FCQPT [3J. In our case F does not 
depend on the number density x of the system as it is 
fixed by condition that the system is situated in QCP 
of FCQPT. Thus, the variational procedure, being ap- 
plied to the functional E[n a (p,T)], gives following form 
for e(p, T) - £ CT (p, T) = e[n a (p, T)] 
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Equations ([3]) and (HJ constitute the closed set for self- 
consistent determination of e a (p,T) and n a (p,T). The 
solution of eq. Q generates the spectrum where the first 
two p-dcrivatives equal zero. Since the first derivative is 
proportional to the reciprocal quasiparticle effective mass 
l/M*, its zero just signifies QCP of FCQPT. The second 
derivative must vanish also. Otherwise e(p) — p, has the 
same sign below and above the Fermi surface, and the 
Landau state becomes unstable [5J[3] • Zeros of these two 
subsequent derivatives mean that the spectrum e(p) has 
an inflection point at pp so that the lowest term of its 
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Taylor expansion is proportional to (p — pf) 3 - In other 
words, close to FCQPT the single - particle spectrum does 
not have customary form vf(j> — Pf), vf is fermi velocity. 

Having solved eqs. ((3]) and (j4]), we substitute their so- 
lution into eq. (H} to obtain field and temperature depen- 
dence of Landau quasiparticlc effective mass. We empha- 
size here, that in our approach the entire temperature and 
magnetic field dependence of the effective mass is brought 
to us by dependencies of e CT (p, T) and n CT (p,T). The sole 
role of Landau amplitude is to bring the system to FCQPT 
point, where Fermi surface alters its topology so that the 
effective mass acquires temperature and field dependence, 
see Ref. and references therein for details. 

Rewriting the quasiparticlc distribution function as 
n a {p,T,B) = n a (p,T = 0, B = 0) + 5n a (p, T, B) yields 
more convenient form for the equation ([T]) 
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Our analysis shows, that near FCQPT the normalized 
solution of eq. ([5]) M^(y = Tn) can be well approximated 
by a simple universal interpolating function. The interpo- 
lation occurs between the LFL (M* cx a + bT 2 ) and NFL 
(M* cx T -2 / 3 ) regimes [2J|5] 
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Here a and b are constants, cq = (1 + c%)/{\ + Ci), c\ and 
C2 are fitting parameters, approximating the Landau am- 
plitude. Note, that our interpolative solution (JB]) is valid 
at low magnetic fields, where spin dependence in Landau 
amplitude and single particle spectrum is not pronounced. 
At high fields, when this dependence is strong and we have 
full subbands spin polarization, this interpolative solution 
is no more valid and we should explicitly solve eq. ([5]) 
with respect to §3§ and (U]). It can be shown that mag- 
netic field B enters Landau equation only in combination 
B^b/T making T max cx B[Lb 0[5]. We conclude that 
under the application of magnetic field the variable 



y = T/T n 
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remains the same and the normalized effective mass is 
again governed by eq. ©• Here B c q is the critical mag- 
netic field driving both HF compound to its magnetic field 
tuned QCP and corresponding Nccl temperature to T = 0. 
In some cases B c q = 0. For example, the HF compound 
CcRu2Si2 has B c q = and shows neither magnetic order- 
ing nor superconductivity [6]. In our simple model B c q is 
taken as a parameter. In what follows, we compute the 
effective mass using eq. © and employ eq. © for qualita- 
tive analysis when considering the system at low magnetic 
fields. 
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Fig. 2: Schematic phase diagram of YbRri2Si2 (© with mag- 
netic field as a control parameter. The vertical and horizontal 
arrows show LFL-NFL and NFL-LFL transitions at fixed B 
and T respectively. At B < B c o the system is in AFM state. 
The width of the transition region T w cx T is shown by the seg- 
ment between two vertical arrows. Inset shows a schematic plot 
of the normalized effective mass versus the normalized temper- 
ature. Transition region, where reaches its maximum at 
T/Tmax = 1, is shown by the hatched area both in the main 
panel and in the inset. 



Now we have everything to construct the schematic 
phase diagram of the HF metal YbRh 2 Si 2 at B <C B* . 
The phase diagram is reported in fig. [5] The magnetic 
field B plays a role of the control parameter, driving the 
system towards its QCP. In our case this QCP is of FC- 
QPT type. The FCQPT peculiarity occurs at B = B c0 , 
yielding new strongly degenerate state at B < B c o. To lift 
this degeneracy, the system forms either superconduct- 
ing (SC) or magnetically ordered (ferromagnetic (FM), 
antifcrromagnetic (AFM) etc) states [2]. In the case of 
YbRh 2 Si2, this state is AFM one pQ. As it follows from 
eqs. (|6]) and (J7J) and seen from fig. [2j at B > B c0 the sys- 
tem is in cither NFL or LFL states. At fixed temperatures 
the increase of B drives the system along the horizontal 
arrow from NFL state to LFL one. On the contrary, at 
fixed magnetic field and raising temperatures the system 
transits along the vertical arrow from LFL state to NFL 
one. The inset to fig. [2J demonstrates the behavior of the 
normalized effective mass = M*/M^ ax versus nor- 
malized temperature y = T/T max following from cq. ©. 
The T -2 / 3 regime is marked as NFL one since (contrary 
to LFL case where the effective mass is constant) the ef- 
fective mass depends strongly on temperature. It is seen 
that temperature region y ~ I signifies a transition regime 
between the LFL behavior with almost constant effective 
mass and NFL one, given by T~ 2 / 3 dependence. Thus, 
temperatures T ~ T max , shown by arrows in the inset and 
main panel, can be regarded as a transition regime be- 
tween LFL and NFL states. It is seen from eq. ([7]) that 
the width of the transition regime T w cx T is proportional 
to (B — B c q). It is shown by the segment between two 
vertical arrows in fig [5J These theoretical results are in 
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good agreement with the experimental facts [UU]. 

Our calculations of the normalized effective mass 
Af^r(TV) at fixed high magnetic field B* are shown by 
the solid line in the main panel of fig. [1] We recollect 
that in this case the quasiparticlcs spins are completely 
polarized. This reveals the above scaling behavior of the 
normalized experimental curves in terms of the normal- 
ized variable y = T/T max (B). It is seen from fig. [T]that 
our calculations deliver a good description of the experi- 
ment pp. Namely, at elevated temperatures (y ~ 1) the 
LFL state first converts into the transition one and then 
disrupts into the NFL state. 

To perceive further the behavior of the system at high 
magnetic fields, in fig. [3] we collect the curves M^(Tiv) 
both at low (symbols in the upper box in fig. [3j and 
high (symbols in the lower box) magnetic fields B. All 
curves have been extracted from the experimental facts 
P~|T5] . It is seen that while at low fields the low-temperature 
ends (Tm ~ 0.1) of the curves completely merge, at high 
fields this is not the case. Moreover, the low-temperature 
asymptotic value of C/T = at low fields is around two 
times more then that at high fields. The physical reason 
for low-field curves merging is that the effective mass does 
not depend on spin variable so that the polarizations of 
subbands with af and ctj, are almost equal to each other. 
This is reflected in our calculations, based on eq. ([BJ for 
low magnetic fields B << B* . The result is shown by the 
dotted line in fig. [3J 

It is also seen from fig. [3J that all low-temperature dif- 
ferences between high- and low field behavior of the nor- 
malized effective mass disappear at high temperatures. In 
other words, while at low temperatures the values of 
for low fields are two times more then those for high fields, 
at temperatures TV > 1 this difference disappear. It is 
seen that these high temperatures lie about the transition 
region, marked by hatched area in the inset to fig. H This 
means that two states (LFL and NFL) separated by the 
transition region are clearly seen in fig. [3J displaying good 
agreement between our calculations (dotted line for low 
fields and thick line at high fields) and the experimental 
points (symbols). 

It is seen from fig. [3J that at high fields B ~ B* , (sym- 
bols in the lower box) the curves M^(Tjy) do not merge in 
the low temperature LFL state. Moreover, their values de- 
crease as B grows representing the full spin polarization of 
the HF band at the highest reached magnetic fields. This 
behavior is opposite to that at low fields. The correspond- 
ing theoretical curve has been generated from the explicit 
numerical solution of eq. ([S]) with respect to eqs. d3J and 
flU . As we have mentioned above, at temperature raising 
all effects of spin polarization smear down, yielding the 
restoration of NFL behavior at T ~ [ibB. Our high-field 
calculations (solid line in fig. [3j reflect the latter fact and 
are also in good agreement with experimental facts. In 
order not to overload fig. [3J with unnecessary details, we 
show the calculations only for the case of the complete 
spin polarization. 
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Fig. 3: Joint behavior of the normalized effective mass at 
low (upper box symbols) and high (lower box symbols) mag- 
netic fields extracted from the specific heat (C/T) measure- 
ments of the YbRh2Si2 [5]. Our low-field calculations are de- 
picted by the dotted line tracing the scaling behavior of MJ. 
Our high-field calculations (solid line) are taken at B ~ B* 
when the quasiparticle band becomes fully polarized. 



Figure H] reports the maxima Mj^ aK (B) of the functions 
in the inset to fig. [T]versus B. The solid line represents our 
approximation for these maxima Af * lax (i3) tx \j \JB — B c q 
calculated within the framework of FCQPT theory [3J[H]. 
It is seen that our calculations are in good agreement with 
the experimental facts in the entire magnetic field domain. 
Such good coincidence indicates that at T ~ ^bB the 
transition regime occurs and the NFL behavior restores 
at high temperatures T ~ 20 K. 
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Fig. 4: The maxima M^^{B) of the functions C/T versus mag- 
netic field B for YbRli2Si2. The points for low [8] (squares) and 
high fields [T] (circles) are shown in the legend. The solid curve 
is approximated by M^ ax (i3) tx d/y/W— B c q, d is a fitting pa- 
rameter. 

In fig. [5l the solid squares and circles denote tempera- 
tures T max at which the maxima of C/T (from the inset 
to fig. [1]) occur. To fit the experimental data [TJ[5] the 
function T max (Ti) = b(B — B c0 ) defined by eq. (JT]) with 
b ~ 0.74 K/T is used. It is seen from fig. [5]that our calcu- 
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lations (solid line) are in accord with experimental facts, 
and we conclude that the transition regime of YbRli2Si2 
is restored at temperatures T ~ hbB. 




Fig. 5: The temperatures T max (_B) at which the maxima of C/T 
in YbRli2Si2 (inset to fig. [TJ are located. Squares correspond 
to low-field case [S] and circles to high fields one Q~|. The solid 
line represents the function T max oc b(B — B c o), b is a fitting 
parameter, see eq. lf?J). 




; Normalized field 

I I I I I 1 — ' ' 

1E-5 1E-4 1E-3 0,01 0,1 

Magnetic field B, arbitrary units 

Fig. 6: The calculated magnetization M(B) (symbols) and 
straight lines which are guides for eye. The intersection of 
the straight lines visualize the kink at the crossover region in 
the fig. [2] The inset: The field dependencies of the normalized 
magnetization M of YbRu 2 Si 2 [TO] at T = 0.05 K. The kink 
(shown by the arrow) is clearly seen at the normalized field 
Bm — B/Bk — 1. The solid curve represents our calculations. 



Consider now the magnetization M(B,T) as a function 
of magnetic field B at fixed temperature T 

M(B,T)= [ x(*,T)dz, (8) 
Jo 

where the magnetic susceptibility \ is given by [3] 

. BM*(B, T) 

X^)- 1 + F a ■ W 

Here, (3 is a constant and F§ is the spin-antisymmetric 
Landau amplitude taken at L = 0. 

Our calculations show that the magnetization exhibits 
a kink at some magnetic field B = B^- The experimen- 
tal magnetization demonstrates the same behavior pTJUTT) . 
We use Bk and M(Bk) to normalize B and M respectively. 
In the normalized variables, there are no coefficients j3 and 
(1 + Fq) so that x M* [2] and we can once more use 
eq. JS]) to calculate the magnetic susceptibility x- The 
normalized magnetization M(B)/M(Bk) both extracted 
from experiment (symbols) and calculated one (solid line) , 
are reported in the inset to fig. El It shows that our cal- 
culations are in good agreement with the experiment. All 
the data exhibit the kink (shown by the arrow) at Bn — 1 
taking place as soon as the system enters the transition re- 
gion. This region corresponds to the magnetic fields where 
the horizontal arrow in fig. [3] crosses the hatched area. To 
illuminate the kink position, in the fig. [6] we present the 
M(B) dependence in logarithmic - logarithmic scale. In 
that case the straight lines show clearly the change of the 
slope (power in logarithmic scale) of M{B) at the kink. 

At magnetic field B ~ B* the quasiparticle band be- 
comes fully polarized and a new kink appears [JJ[T2] . We 



call this kink as the second one. Our calculations of 
the normalized magnetization (line) and the experimen- 
tal points (squares) are shown in fig. [71 In that case both 
the magnetization and the field are normalized by the cor- 
responding values at the second kink position. 
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Fig. 7: The normalized magnetization M nor as a function of 
the normalized magnetic field Bn- Line represents our cal- 
culations and squares represent experimental points [JJ. The 
linear dependence M nor (Bjv) is marked by the arrows. The 
inset demonstrates the experimental data (stars) with the sub- 
tracted high-field linear part. Our calculations is shown by the 
solid line. 

In the fig. [JJ we plot our theoretical normalized (in the 
second kink point) magnetization along with experimen- 
tal one. Good coincidence is seen everywhere except the 
high- field part at Bn > 1. Here, the experimental nor- 
malized magnetization M nor exhibits a linear dependence 
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on Bn (marked by two arrows), while the calculated mag- [15] Rourke P.M.C. et al, Phys. Rev. Lett., 101 (2008) 
netization is approximately constant. Such a behavior is 237205. 
the intrinsic shortcoming of the HF liquid model that ac- 
counts for only heavy electrons and omits the conduction 
electrons of other kind [T51ll4j . Thus, we can consider the 
high-field (at Bm > 1) part of the magnetization as the 
contribution which is not included in our theory. To sep- 
arate this contribution from the experimental magnetiza- 
tion curve, we (numerically) differentiate it, then subtract 
constant part at Bn > 1 and integrate back the resulting 
curve. The coincidence between our calculations depicted 
by the solid curve and processed experimental data shown 
by the stars is reported in the inset to fig. [JJ As we can see 
now, the coincidence between the theory and experiment is 
good in the entire magnetic field domain. Taking into ac- 
count the obtained results displayed in figs. [3l []] O [6] and 
we conclude that the HF system of YbRu2Si2 evolves 
continuously under the application of magnetic field. This 
fact is in agreement with experimental observations [15] . 

To summarize, here we have analyzed the thermody- 
namic properties of YbRli2Si2 at both low and high mag- 
netic fields. Our calculations allow us to conclude that in 
magnetic field the HF system of YbRli2Si2 evolves con- 
tinuously without a metamagnetic transition and possi- 
ble localization of heavy 4/ electrons. Under the applica- 
tion of magnetic field at low temperatures, the HF system 
demonstrates the LFL behavior, while at elevated temper- 
atures the system enters the transition region followed by 
the NFL behavior. Our calculations are in good agree- 
ment with experimental facts in the entire temperature 
and magnetic field domains under consideration. 

This work was supported in part by the RFBR grant 
No. 09-02-00056. 
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